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A NEW ISOSCELES-TRIANGLE SOLUTION OF THE THREE 

BODY PROBLEM* 

BY 

DANIEL BUCHANAN 

1. Introduction 

The known isosceles-triangle solutions of the problem of three bodies are 
the periodic solutions in which two of the bodies have equal finite masses, 
while the third body moves in a straight line and remains equidistant from 
the other bodies. In my Chicago dissertation! three cases of the isosceles- 
triangle solutions were discussed. In the first case the finite bodies move 
in a circle and the third body is infinitesimal. The solution of this case was 
first obtained by PavaniniJ by means of elliptic integrals, and was later de- 
veloped independently by MacMillan.§ The solution which was obtained 
in my dissertation is analogous to that obtained by MacMillan. In the 
second case, not treated by Pavanini and MacMillan, the finite bodies move 
in ellipses and the third body is infinitesimal. In the third case, all three 
bodies are finite. The period of the motion in case III is the same as that 
in case I, and the solution reduces to the solution of the first case when the 
third body becomes infinitesimal. 

The present paper is an extension of the second case of my dissertation to 
that in which all three masses are finite. Periodic orbits are here shown to 
exist, and a practical method is given of computing them. These solutions 
have the same period as those of case II of my dissertation, and reduce to them 
when the third body becomes infinitesimal. 

2. The differential equations of motion 

Let mi and mi represent the two finite bodies of mass m , and let the third 
body be denoted by /x . The origin of coordinates will be taken at the center 
of mass of the system. The plane perpendicular to the initial motion of n 

* Presented to the Society (Chicago), December 28, 1914. 
t Moulton's Periodic Orbits, Chap. X. 

t A n n a 1 i d i M a t e m a t i c a , ser. 3, vol. 13 (1907), pp. 179-202. 
§Astronomical Journal, Nos. 625-626 (1911). 
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will be taken as the ^-plane. Then let the coordinates of mi, mi, and ju be 
denoted by £1, i\x, ft; £ 2 , 172, ft; and £, 77, f respectively. 
When ju = the mean angular motion of the finite bodies is* 

k V2m 
(!) n = ~3"/2- > 

where k is the gravitational constant and a is the mean distance from mi to mi . 
Now let the unit of mass be so chosen that m = 1/2 , and let the linear and 
time units be so determined that k and a are each unity. With the units 
thus chosen n = 1 . 

We shall suppose that ix is started from the origin of coordinates at the 
initial time to when mi and mi are at apsides of their orbits, which can be 
assumed, without loss of geometric generality, to lie on the £-axis. When 
expressed analytically, these initial conditions are 

£ = i? = £ = £' = 17' = i?i = i?2 = ft = ft = at t = t , 

fc(*o)--€i(*o)=f(l-«), *i«o)=-£(*o), i?'i«o)=-i?i«o), 

where e is the eccentricity of the ellipses in which mi and mi move when 
H = . Under these initial conditions the orbits described by mi and m 2 will 
be symmetrical; that is, 

(2) £1 = - £ 2 , 171 =-172, ft = ft. 

On making use of (2) in the center of gravity equations, which are 

mi £1 + m 2 £2 + m£ = 0, mi rn + m 2 172 + A"7 = 0, 

mi ft + m 2 ft + Mf = 0, 
we obtain 

(3) £-17 = 0, ft + Mr = o. 

Hence the third body always remains on the ftaxis. 

The differential equations of motion when simplified by (3) become 

t" - ^ vh 

£> — 0»3 



8r 3 [r 2 + (1 + M) 2 r 2 ] 3/2 ' 



W 171 - gr 3 [ r 2 + (1 +M )2 f 2]3/2> 

-2m(l+/i)r 
f [r 2 + (1 + M ) 2 r 2 ] 3/2 ' 

where r 2 = £1 + r\\ . When these equations have been integrated, the solu- 
* Moulton's Celestial Mechanics, p. 142. 
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tions for £ 2 , tyz> fi» and £"2 can be obtained from (2) and (3). The numerical 
value of m = 1/2 has been substituted in the first two equations of (4) but 
not in the last equation. The reason for this will appear in the next section. 
On transforming to polar coordinates by the substitutions 

£1 = r cos v , 571 = r sin v , f = z , 

the differential equations become 

„M2 ■ 1 » r 



r"-r(v'y + ^ 2 = 



8r> [r 2 + (1 + m) 2 2 2 ] 3/2 ' 

(5) rv" + 2rV = 0, 

„ _ 2m ( 1 + M ) z 

2 " [f«+ (1 + M) 2 2 2 ] 3/2 ' 

3. The equations of variation 

For ix = 0, the first two equations of (5) reduce to those of the two-body 
problem, and their solutions are* 



r 



= r = K 1 — e cos n(t — t ) + » {1 — cos 2n(t — t )}+ • •• , 



( 6 ) 5 
v = v = n{t — t ) + 2e sin n(t — t ) + je 2 sin 2n(t — t ) + 

In these equations both a and n are unity, while t can be chosen to be zero 
without loss of generality, but these numerical values are not substituted at 
present as we wish later to differentiate r with respect to t . 
The second equation of (5) admits the integral 

(7) r 2 v' = c. 
When n = the third equation of (5) becomes 

2mz 



(8) z" = - 



[r 2 -|-z 2 ] 3/2 ' 



where r has the value r obtained in (6). Now instead of putting m = 1/2 
in the numerator of (8), as was done in the first two equations of (4), we make 
the substitution 

m = mo + \, 

where mo is considered as a constant and X as a variable parameter. For 
the purpose of computation we must substitute X = m — mo in the final 
results. 



* Moulton'8 Celestial Mechanics, pp. 154, 155. 
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It was shown in case II of my dissertation that periodic solutions of equa- 
tion (8) exist, and are power series in X 1 / 2 of the form 

zo = ± X 1/2 — p= sin an {t — U) + power series in e + 
L 3 Vm J 

where each z 2 /+i is a power series in e of the form 

22/+1 = Z|Z £<##.« sin n{(2* + l)cr ± i}(t - t )\e*. 
K =o L 1=0 <=o J 

The cttj, « are rea l constants and a is a power series in e of the form 

3A/m (l -16otq) 2 ^ 
<7 = 4Vm + ^^ e 2 +.... 

As power series similar to zo in (9) occur frequently in the sequel, we shall 
call them ^^-sine-series. 

Since a is a continuous function of e, values of e can be chosen so that a 
converges to a rational fraction N/v, where N and v are integers relatively 
prime. The period of the solutions (9) is then 

do) r = 2,x = 2 -^. 

In § 171 of my dissertation,* a reference was given to show that the two 
solutions in (9) represent the same orbit if v is even, but in the one the in- 
finitesimal body is half a period ahead of its position in the other. If v is 
odd the orbits for plus and minus are geometrically distinct. In the sequel 
we shall take only the plus sign in (9) as the solutions obtained from the minus 
sign are entirely analogous to those obtained when the plus sign is taken. 

In order to find the solutions of (5) when \x is not zero, we make the sub- 
stitutions 

(11) r = r + p, c = c + ^M. z = z + w, 

where p and w vanish with ju, and d is an arbitrary constant. From (7) it 
follows that 

Co = rl v' . 

When v' is eliminated from the first equation of (5) by the integral (7), and 
the substitutions (11) are made, the terms in the first equation of (5) inde- 
pendent of n are zero. Hence 

(12) d =r' 'rl + ±r . 



Moulton's Periodic Orbits. 
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The differential equations determining p and w then become 

p" + P(e,t)p = Po + Piix + P 2 ix 2 + ■■-, 
(13) 

w" + W(\,e,i)w = W + W lt x + W if i 2 + ••-, 

where P , W , Pj, Wj (j = 0, 1 , 2, • • • ) have the following properties: 

(1) P is a power series in e whose coefficients are sums of cosines of multi- 
ples of t . In so far as the computations have been made, we have 

P(e,t) = 1 + 6ecos< + e 2 (3 + 12 cos 20 + •••• 

(2) W is a power series in X whose coefficients are power series in e . The 
coefficients of the various powers of e have the form 



Vc fc ,i cos (2Jc<t ± l)t, 



where Jc and I are positive integers or zero, and Ck,i are real constants. As 
power series of this kind occur frequently in the sequel we shall call them X- 
cosine-series. 

(3) Po is a power series in p beginning with the term in p 2 . The coefficients 
of the various powers of p are power series in e similar in form to P . 

(4) Pj (j = 1,2,3, • • • ) are power series in p and w. The coefficients 
of the terms in which w occurs in even degrees are X-cosine-series. The 
coefficients of the terms in which w occurs in odd degrees are X 1/2 -sine-series. 

(5) Wj (j = 0,1,2, •••) are power series in p and w . The coefficients of 
the various powers of p and w are X 1/2 -sine-series if w occurs in even degrees, 
and X-cosine-series if w occurs in odd degrees. 

If we put the left members in (13) equal to zero, we obtain the equations 
of variation, which are 

(14) (a) p" + P{e,t)p = 0, (6) w" + W (\, e, t)w = 0. 

4. The solutions of the equations of variation 

The generating solutions of the equations of variation are ro and so in (6) 
and (9) respectively. Each of these two generating solutions contains the 
arbitrary constant t , and hence one set of solutions of (14) is obtained by 
taking the first partial derivatives of ro and 2o with respect to to* One 
solution of the first equation of (14) is therefore 

(15) ^= -^[Ant + eAn2t+ ■■■] = --^(e,t), 

and as we multiply this solution later by an arbitrary constant, we may dis- 
regard the factor — e/2 . 

* Poincar6's Les Methodes NouvelUs de la Micanique Cileste, vol. 1, Chap. IV. Moulton's 
Periodic Orbits, §§ 32, 33. 

Trans. Amer. Math. Soc. 18 
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When the generating solution contains an arbitrary constant one of the 
characteristic exponents of the solutions of the equations of variation is zero, 
and since the sum of the characteristic exponents is zero both these exponents 
are zero. In this case one of the solutions is a periodic function while the 
other is the sum of a periodic function and t times the former solution. That 
is, the second solution of (14a) is 

where 8 and A are power series in e , the former having periodic coefficients, 
the latter having constant coefficients. In order to determine 8 and A we 
substitute p = 8 -f- Afy in (14a) and integrate as power series in e . We 
choose the initial values 0(O)=1,0'(O) = O and determine A so that 8 
shall be periodic. The computation gives 

8 = cos t — e ( 3 — 2 cos t — cos 2t ) 

- e 2 (6 - V- cos t - 2 cos 2t - f cos 30 + • • • , 

A = ye 2 = [ 3 + power series in e ] e 2 . 

Hence the two solutions of (14a) are $ and 8 + ye 2 fy . The determinant 
of these two solutions is a constant,* and its value at t = is 

(16) D= -8(0)V(0) = -(l+2e+ •••)• 

This determinant is different from zero for e = and therefore remains 
different from zero for | e | sufficiently small. Consequently the two solutions 
(15) constitute a fundamental set and the general solution of the first equation 
of (14) is 

(17) p = A 1 ^ + A 2 [8 + ye 2 t^], 

where Ai and A2 are the constants of integration. The functions ^ and 8 
are periodic with the period 2ir and a posteriori with the period T . 
Similarly, one of the solutions of the second equation of (14) is 

(18) ^ = V^(X,«,(), 

where is a power series in X similar to a X-cosine-series except that it has 
odd multiples of cr in the arguments of the cosines instead of even multiples* 
By computation it is found that 

(19) 0(0)=-^+..., 0'(O)=O. 

3-Vmo 

The second solution of (146) could be found in the same way as the second 



' Moulton's Periodic Orbits, § 18. 
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solution of (14a), but we shall determine it by a shorter method. We may 

consider X as an arbitrary constant and in this case the other solution of (146) 

is 

„ m dz /dz \ dzpdn 

{M) d\ -\d\) + dnd\' 

where (dz /dX) denotes that the differentiation is performed only in so far 
as X enters zo explicitly. Hence 



(dzo\ _ dz d(X 1/2 ) 1 

\d\) d (X 1 ' 2 



!/ 2 ) 3X X 1/2X ' 

where x is a power series in X with coefficients which are similar in form to 
the coefficients in a X 1/2 -sine-series. The values of x (0) and x' (0) are found 
by computation to be 

f V2<r 1 

(21) x (0) = 0, x' (0) = — ~ + power series in X . 

L6A/m J 

The partial derivative dzo/dn has the value 

and since 

k^2m fcV2(m + X) 
n ~ a 3 ' 2 ~~ a 3 ' 2 

then dn/d\ = 1 . Therefore the second solution of (146) is 

-^ITilx +X<0]. 
The two solutions of (146) then take the form 

(22) X 1/2 and -^-Jx + X^], 

where <f> and x are periodic with the period T . As in the solution of (14a) 
we can omit the constants X 1/2 and — X -1/2 , and consider $ and x + X£$ as 
the two solutions. The determinant of these solutions is a constant,* and 
at t = the value is 

(23) A -0(O)[ X '(O)+X^(O)] = — +••-, 

a power series in X in which the coefficients are power series in e . The deter- 



' Moulton's Periodic Orbits, § 18. 
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minant A is different from zero for e = X = , and therefore remains different 
from zero for \e\ and |\| sufficiently small. Consequently the two solutions 
<f> and x + X20 constitute a fundamental set and the general solution of (146) is 

(24) w = Bi* + Atx+W0], 
where B\ and B 2 are the constants of integration. 

5. Existence of symmetrical peeiodic orbits 
Let us take the initial conditions 

(25) p(0) = «i, 2>'(0) = 0, w(0) =0, w'(0)-« t . 

With these initial conditions it can be shown from the properties of the differ- 
ential equations (13) that p is even in t and w is odd in t . Hence sufficient 
conditions that p and w shall be periodic with the period T are 



(26) 



P'(f) = 0, «(f)-0. 



In order to establish the existence of symmetrical periodic orbits, we integrate 
(13) as power series in a\, az, and ju, subject to the initial conditions (25), 
and impose the periodicity conditions (26) upon the solutions. It will be 
shown that (26) can be satisfied by a determination of «i and a 2 as power 
series in n . 

The differential equations from which the linear terms in «i and a 2 are 
obtained are the same as the equations of variation. The solutions of these 
equations which satisfy the initial conditions (25) are 

«i tn , ,.,, «2[x + M<t>] 

p = eW) le + yem ' w = ITW+MWY 

When the terms in n and higher degree terms in n , «i , and a 2 are taken in 
(13), the solutions satisfying (25) have the form 

+ terms in y. and higher degree terms in n , «i , and «2 , 

(27) « 2 [x+A^] 

w = 



x'(0)+x<M0)] 

+ terms in /* and higher degree terms in ju , «i , and a 2 . 



0= ai 



(28) 

«2 
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After imposing the conditions (26) upon (27), we obtain 

0(o)KD + ^'(f)+Wf)] 

+ terms in ju and higher degree terms in ju , «i , and «2 , 

[»(f)+*HI)l 

X ' + X0(O) 

+ terms in ju and higher degree terms in n , «i , and a% . 
From the form of \f/ , 6 and x it follows that 

♦(!)-*(!)-*(?)-■ 

and therefore the determinant of the linear terms in a\ and a 2 in (28) 

T^rvfj) 0(f)^(O)0(O) 



(29) 5 = - 



4Z>A 



The value of 5 depends upon the parity of v and N in T = 2vir = 22W/<r. 
The determinants D and A are both different from zero for \e\ and |X| suf- 
ficiently small and are both finite. Hence 8 can become zero only when the 
numerator vanishes. We shall show that the numerator of (29) is different 
from zero provided that e and X are not zero but sufficiently small numerically. 
Suppose v is even, then N is odd, since v and N are relatively prime, and from 
the form of \f/ and 4> it follows that 

*'(f) = *'(0)*0, *(f)=-*(0)*0. 
If v is odd then 

*'(f) = ~ l + 2e +'--> 



and 



fT\ V2 . V2 

<j> I — J = — = + e [power series in e] or — — -=■ + e [power series in e] , 
\ ■* / 3 Vrao 3 Vm 

according as iV is even or odd. Therefore 8 is different from zero whatever 
be the parity of v and N provided \e\ and |X| are sufficiently small. Conse- 
quently equations (28) can be solved for «i and ai as power series in y. which 
vanish with n and converge for |ju| sufficiently small. Since the initial 
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values of p and w are power series in n , the solutions will have the form 

00 00 

(30) p = 2 Pi m\ «> = H w% M f , 

where each p»- and Wi is periodic with the period T . The constant d which 
appears first in (11) enters these solutions as an arbitrary, but when a value 
is assigned to d the solutions in (30) are unique. 

6. Direct construction of the symmetrical periodic orbits 

Let us substitute (30) in (13) and equate the coefficients of the various 
powers of /* . We thus obtain a series of differential equations which can be 
integrated step by step. The constants of integration appearing at each 
stage are to be determined so that the solutions shall be periodic and satisfy 
the symmetrical initial conditions 

p'(0) = 0, w(0) = 0. 

When these symmetrical initial conditions are imposed on (30) we obtain 

(3i) p;(o) = o, w*(o) = o (t = i, •••«). 

The differential equations for the terms in ju are 

(32) (a) p[' + P(e,t)p 1 = PV > (b) w'{ + W (X, e, t) Wl = W™ . 

The left members of these equations are the same as the equations of vari- 
ation. The right member P (1) is a X-cosine-series and is independent of Wi . 
The right member of (326) depends on pi. The terms of W a) which are 
independent of pi form a X 1/2 -sine-series and the remaining terms are also a 
X 1 /2 -sine-series multiplied by p\ . 

The complementary function of (32a) is 

(33) pi = a\ l) + + ayie + yeHt], 

where a^ and a ( 2 1} are the constants of integration. To obtain the particular 

integral we employ the method of the variation of parameters. By this 

method we have 

(a[»y+ + (dPy[6 + yeHt] = 0, 
(34) 

«)' V + (WY [6 + ye* W + *) ] = P0> . 

The determinant of the coefficients of (a^)' and (a^)' is D, the same as 
in (16), and is different from zero. Hence 

(35) w- ' e+ t" f '° - w-*£- 
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When these equations are integrated and the results are substituted in (33), 
we obtain the solution 

(36) p 1 = A[ 1 ^+A ( p[d+ye 2 m+C^-t^[L^ + \ v ' 2 e x L[ l) or \"e 2N L^], 

according as v is even or odd respectively. The terms in (36) not already 
defined have the following properties: A^ and A ( P are the constants of 
integration; C (1) is a X-cosine-series; Z^ 1 ' and Z4 1 ' are power series in X with 
coefficients which are constant power series in e, and i'o 1 ' is a power series 
in e with constant coefficients. In order that pi shall satisfy (31) and shall be 
periodic, the constants of integration must have the values 

(37) AV-0, AP=— 2 lX ) +-[\ l " 2 e N - 2 LP or \" e 2N ~ 2 W ] , 

ye y 

according as v is even or odd. Hence the solution of (32a) which is periodic 
and satisfies (31) is 

(38) Pi-4Ci(0, 

where C\ ( t) is a X-cosine-series. 

With pi determined as in (38), the right member of (326) takes the form of 
a X 1 /2 -sine-series multiplied by 1/e 2 . The complete solution of (326) is obtained 
in the same way as (36) was obtained. The complementary function of 
(326) is 

(39) w>i = 6 ( 1 1) <£ + 6 ( 2 1) [ x + X*0], 

and the equations analogous to (34) are 

(6 ( 1 1 0'^ + (& ( 2 1) )'[x + X^] = 0, 

(#»)' 0' + m'[ x ' + X(*0' + 0)] = W». 

The determinant of the coefficients of (6 ( 1 1) )' and (6 ( 2 1) )' is A, the same as 
in (23), and is different from zero. Hence 

(40) <„,,.■. faLfaJMig , W y.Hp. 

When these equations are integrated and the results are substituted in (39), 
the general solution of (326) is found to be 

(41) Wl = M l> <t> + W \ x + w* ] + f 2 s« - ^ JJf a) t* , 

where B^ and B^ are the constants of integration; S (1) is a X^-sine-series, 
and M w has the same form as Z4 1} or Z4 1 ' in (36). In order that Wi shall 
satisfy the initial conditions (31) and be periodic, the constants of integration 
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must have the values 

(42) 2^ = 0, Bf = ^. 
Hence the desired solution of (326) is 

(43) «i=-^Si(0, 

where Si ( t) is a X 1 / 2 -sine-series. 

The remaining steps of the integration can be carried on in the same way, 
and the solutions are almost entirely similar to those obtained at the first step. 
The presence of the terms in 1/e 2 and l/e 2 X in (38) and (43) introduces some 
slight changes in the forms of the solutions, and it is necessary to consider 
the differential equations up to the terms in ju 3 in order to discover the law 
according to which the remaining solutions proceed. 

The terms in ju 2 yield the differential equations 

(44) (a) pi +P(e, t)p 2 = P (2) , (6) w 2 + W(\, e, t)w 2 = W*> . 

The right member of (44a) is a X-cosine-series multiplied by 1/e, and the 
complete solution of the equation is 

(45) p 2 = A? f + Af [d + yeHH+i C« - 4 Z« f* , 

where A ( ? and A^ are the constants of integration; C (2) is a X-cosine-series; 
and i (2) is a power series in X with coefficients which are constant power series 
in e. The numerical values of the coefficients in i (2) vary with the parity 
of v , but the series has the same form whether v is even or odd. In order 
that pi shall be periodic and satisfy the conditions (31) the constants of 
integration must have the values 

(46) A?> = 0, Af = ^, 
and the desired solution for p 2 then takes the form 

Pi =- 4 C 2 (t), 

where Ci ( t) is a X-cosine-series. 

The right member of (446) is completely known when p 2 has been deter- 
mined and is a X^ 2 -sine-series multiplied by 1/e 4 X . The complete solution 
of (446) is readily found to be 

(47) w 2 = Bf + Bf [ X + W0 ] + ^ « (2) - ^I7i M& Uj> , 

where the undefined terms in (47) denote functions similar to those represented 
by the corresponding terms in (41). If the initial and periodicity conditions 
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are to be satisfied then 

Hence the desired solution for w 2 is 

where S 2 (t) is a X 1/2 -sine-series. 
The differential equations obtained from the terms in n s are 

(48) (a) pi +P(e,t)p z = P (3 \ (b) w'i + W (X, e, t)w 3 = W™ . 

The right member P (3) is a X-cosine-series multiplied by l/e 4 X, and the com- 
plete solution of (48a) is 

(49) p 3 = A? + + Af [e + yeHxf,] + e lr C< 3 > - -^ Z< 3 > ty> , 

6 A € A 

where the undefined terms denote functions similar to those represented by 
the corresponding terms in (45). If the initial and periodicity conditions are 
to be satisfied, then 

J^ 
ye e \ J 
and the desired solution for p 3 becomes 



4?> = 0, Af>'zk LW > 



where C 3 ( t) is a X-cosine-series. 

From the forms of pz and w 2 it follows that W (s) is a X 1/2 -sine-series multi- 
plied by 1/e 6 X 2 . Hence the solution for w 3 will be a X^-sine-series multiplied 
by 1/e 6 X 3 . By an induction to the general term it can readily be shown that 
Pj and Wj have the forms 

Pj=fTsJ=,Ci(t) U-2,3,— -), 

Wj = -2JxSj(t) = 1,2,... •), 

where each Cy (t) is a X-cosine-series and each Sj (t) is a X^-sine-series. 

7. PEOOF THAT ALL PERIODIC ORBITS REPRESENTED BY POWER SERIES IN fl 

ARE SYMMETRICAL 

The usual method of proving that all periodic orbits are symmetrical is to 
assume the initial conditions 

(49) p(0) = oi, p'(0)~a 2> w(0)-o», w'(0) = 0, 
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and integrate equations (13) as power series in a, (i = 1, 2, 3) and ju. Suf- 
ficient conditions that the solutions shall be periodic are 

p(T)-p(0) =0, p'{T) -p'(0) =0, w{T) -w(0) = 0, 
(50) 

w'(r) - w'(0) = 0. 

If there exists an integral of the differential equations (5) it could be shown 
by means of this integral that the last condition of (50) would be reduntant. 
Then the three independent conditions of (50) could be satisfied by a proper 
determination of a,- ( i = 1,2,3) as power series in ju provided a certain 
determinant did not vanish. Thus periodic solutions would be shown to 
exist and they would be unique for any assigned value of d which would enter 
the solutions through (11). Since the symmetrical solutions are also unique 
under the same circumstance and since the general orbits having the initial 
conditions (49) would include those which are symmetrical, all the solutions 
would therefore be symmetrical. 

In the problem under consideration no integral of equations (5) is known 
to exist and we can not use the foregoing method. We show, however, 
from the construction of the solutions that the symmetrical initial conditions 
are a consequence of the periodicity conditions and therefore that all the 
orbits represented by solutions which are power series in n with periodic 
coefficients are symmetrical. 

If the symmetrical initial conditions (31) are not imposed upon (36), the 
periodic solution for pi is 
(51) p 1 = AV4> + C 1 (t). 

When this value of pi is substituted in the right member of the differential 
equations defining Wi we obtain 

w" + W{\, e, t) Wl = W™ + A?' W®, 

where W m is similar to a X 1/2 -sine-series except that it has cosines instead of 
sines. The complete solution of this equation is the same as (41) with the 
additional terms 

^'i 1 ' [ periodic terms + terms in f 4> and t\ ] • 

The only way in which this solution can be made periodic is by putting A™ = 
and then determining B^ as in (42). Hence the periodic solution for wi 
which does not satisfy the condition wi ( ) = becomes 

Wl = 5 ( 1 1) 0+^Si(O, 
where B ( P remains undetermined at this step. In order to obtain the value 
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of Bp we must consider the differential equation which defines p 2 . It has 

the form 

(52) pi' + P ( e , t) p 2 = P» + Bf} ) P (2) , 

where P (2) is similar to a X 1/2 -sine-series except that it has even multiples of a 
instead of odd. The complete solution of (52) is the same as (45) with the 
additional terms 

B'P [ periodic terms + terms in ft if/ and tO ] . 

The only way in which p 2 can be made periodic is by putting Bf-^ = and then 
determining A^ in the same way as in (46). 

Thus the integrations of the differential equations arising at each step can 
be obtained in the same way and the process continued as far as desired. 
The constant of integration which is not determined by the condition that pj 
shall be periodic must be zero in order that wj shall be periodic; and similarly 
the constant of integration which is not determined by the condition that wj 
shall be periodic must be zero in order that pj+i shall be periodic. The solutions 
thus obtained are the same as those obtained for the symmetrical orbits. 
Hence the symmetrical initial conditions (31) are a consequence of the peri- 
odicity conditions and all periodic orbits represented by power series in ju are 
symmetrical. 

8. Integration of equation (7) 

There remains one differential equation which has not been integrated, 
viz., 
(7) r 2 »' = c. 

In order to integrate this equation we make in (7) the substitutions 

v = v + q, r = ro + p, c = c a + dn, 

where q vanishes with n , and p has the solution obtained in § 6. Then equa- 
tion (7) becomes 
(54) q' = QiH + Q2H 2 + ■■■ +Qjl* i + ■■-, 

where each Q 3 is similar in form to pj, that is, each Qj is a X-cosine-series 
multiplied by l/e 2y X y ~ 2 . As the right side of (54) contains constant terms, 
the integration will yield terms in t unless these terms can be made to vanish. 
Now the constant d appears in these terms but as the constant part of the 
coefficient of n does not vanish when d = 0, it is impossible to determine d 
so as to eliminate the constant terms in (54). Hence the integration of (54) 
gives 

00 

q = YjUj + %]M y + go, 
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where each qj is similar in form to Qj except that it has sines instead of 

cosines, each qj is similar to Qj(0), and go is the constant of integration. 

Since the equal masses are started from the £-axis at the initial time then 

v ( ) = . Since v ( ) = it follows that q ( ) = and as each qj ( ) = 

then go = . 

This completes the integration of the differential equations defining the 

motion of the three finite bodies mi , m 2 , and n . The constant of integration 

arising from the integral of areas remains undetermined, but the remaining 

constants of integration can be determined by the periodicity conditions 

alone. The mass ju is restricted in magnitude because of certain convergence 

conditions but it can be increased step by step by making the analytical 

continuation of the solutions already obtained, provided the series do not 

pass through any singularities in the intervals. This can be done by the 

method already developed. 

Queen's University, Kingston, Ontario, 
August 17, 1914. 



